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Abstract 

We consider a variation of 0(A^)-symmetric vector models in which the vector com- 
ponents are Grassmann numbers. We show that these theories generate the same sort 
of random polymer models as the 0{N) vector models and that they lie in the same 
universality class in the large- limit. We explicitly construct the double-scaling limit 
of the theory and show that the genus expansion is an alternating Borel summable series 
that otherwise coincides with the topological expansion of the bosonic models. We also 
show how the fermionic nature of these models leads to an explicit solution even at 
finite- for the generating functions of the number of random polymer configurations. 
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The statistical mechanics of randomly branching polymers (sometimes called discrete fil- 
amentary surfaces) is of interest in condensed matter physics where the connected polymer 
chains can be thought of as describing molecules. Random polymers have also been of inter- 
est as examples of random geometry systems which in many cases are exactly solvable and 
which have been argued to represent a certain dimensionally-reduced phase of Polyakov string 
theory in target space dimensions D > 1 (^0. They therefore serve as toy models for more 
complicated higher dimensional problems such as statistical models of discretized surfaces 
which have been studied in the context of string theory and lower dimensional quantum grav- 
ity (see for a review). 0(A^)-symmetric vector field theories provide non-perturbative 
models of randomly branching chains where the order in the expansion coincides with 
the genus, or number of loops in the molecules of the ensemble. These models exhibit phase 
transitions in the large- limit at which infinitely long polymers dominate the statistical sum 
and a continuum limit analogous to that found in the matrix model representation of random 
surface theories is reached The simplicity of the random polymers as compared to 

random surfaces allows a more explicit solution, even in dimensions D > 1, where some of the 
ideas about scaling and other critical behaviour can be tested. These simpler polymer struc- 
tures are reflected in the linearity of the vector field theories in contrast to the non-linearity 
of matrix ones. The multicritical series generated by these models give generalized statistical 
systems which interpolate between the Cayley tree at one end and the ordinary random walk 
at the other Some supersymmetric generalizations of the 0{N) vector model have been 
studied in [|10|] . 

In this Letter we shall show that a vector theory with purely fermionic degrees of free- 
dom can also be used to represent random polymers. This model exhibits the same critical 
behaviour in the large- limit as the 0{N) vector model, but it has a structure which is in 
some respects simpler. In particular, we will show that its genus expansion is an alternating 
series which is Borel summable and that it represents a rare example of a random geometry 
theory whose explicit solution can be written down even at finite A^. It therefore provides an 
explicit, well-defined generating function for a given number of polymer configurations in a 
random polymer system. It can also be combined with other polymer models to generate new 
types of generating functions for random surface theories. For example, it can be combined 
with the 0{N) vector model to study random polymer theories with either only even or only 
odd genera. 

Let us start by reviewing some features of the 0{N) vector model with partition function 



Zs{t,g-N) = 




where the integration is over IR^. The model (|l]) is invariant under the orthogonal transfor- 
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mation 0j — > ^f=i ^ij't'j^ S E 0{N). This symmetry restricts the observables of the theory to 
those which are functions of 0^ = J^i'P'i- When the couphng constants t and g are positive, 
it is straightforward to show that the formal expression (|I]) counts the number of randomly 
branching polymers, both those with a tree-like structure and those with arbitrarily many 
loops. The fact that the integral is divergent is a reflection of the divergence of the statistical 
sum. 

The statistical sum over polymers coincides with the expansion of the free energy 
in Feynman diagrams. The propagator is 
the vertex is 

{4>i4>j4>k4>i) s = j2t^ i^ij^ki + Sii6jk + SikSji) (4) 

and the free energy is the sum of all connected diagrams with four-point couplings (Figs. 1 and 
2). The dual graphs to these Feynman diagrams, defined by associating a vertex (molecule) 
in the center of each of the scalar loops and lines (bonds) connecting vertices by crossing 
each of the Feynman 4-point couplings, are the random walk diagrams shown in Fig. 2. The 
number of molecules n appears as the power in A^"~^ while the number of bonds b is given 
by the power in {g/Nt^)^ |^, ^. This identifies = with the fugacity fi of the polymer 
and ^ = e~^ with the length L of the branched chain (the usual action terms for a random 
walk model) The vector model partition function is therefore the generating function 
for the number of polymer configurations with b bonds and i = b — n + 1 loops. Note that 
this sum includes the self-bonding polymers which are generated by Wick contracting several 
propagators into single loops (as opposed to multi-loops) and occur in the expansion only for 
i > 1 (Fig. 2). 



i J 

= ^ij ¥ = iF i5ij5ki + 5ii5jk + SikSji) 




Figure 1: Feynman rules for the 0{N) vector model. 




Figure 2: Feynman diagram expansion of the D = quantum field theory (|^). The thick hues 
and vertices represent the Feynman graphs, while the thin lines and vertices represent the 
dual graphs which form polymer networks (or random walks). The first set of diagrams 
represent those types of branching polymers associated with the higher orders of in vertex 
contractions such as {{4>'^)'^)s = {^"^ + 2N)/4:t'^ in Fig. 1 (so that n > 1 and bonds connect 
different molecules), while the second set of diagrams representing the self-bonding polymers 
are associated with the lower powers of N (so that i = b). 



From an analytic point of view, the perturbative expansion of the partition function 

Zsit,g;N) ^ 1 f 9 Y / f ^2\^n\ 



is completely determined by the Gaussian moments 



^ IJhdMfr^ _ .N/2 ( ^\\-N/2 



s 



as 

Na 



Zs{t, 9;N) _^ {N + 4n-2)\\ f j_ . 



Here we have introduced an "ultraviolet" cutoff A^a G to make the partition function 
well-defined. In the limit A^a oo, the series (^ is a non-Borel summable asymptotic series 
reflecting the divergence of the original integral and also the divergence of the statistical sum. 
As is familiar from the study of random surfaces 0], even though the series is divergent, if 
arranged as a power series in rather than g, the terms in this series are individually 

convergent and it is the sum over genera £ which is asymptotic With the cutoff 

A^A in (^, the partition function is an analytic function of A^, but it is only well-defined at 
N = oo when this ultraviolet cutoff is removed. 

The large A^ expansion is a saddle point computation of the integral ([1|). After integra- 
tion over angular variables, the partition function (|I]) can be written in terms of the radial 
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coordinate in Euclidean A^-space as 

In the infinite- iV limit, which counts the tree-graphs {£ = 0), the integral (||) can be evaluated 
using the saddle-point approximation. Rescaling (f)/\/N the stationary condition for the 
effective action -Nt(j)^ + gNcf)'^ + Nlogcf) in (D is 

The solution of which is regular at = and which minimizes the effective action is 




(10) 

The tree-level free energy is then the value of the effective action in (B) evaluated at the 



saddle-point (|T0[) 



.(0) - 1 . „ Ill t / rz -\ . / t 1 



(11) 

This free energy becomes non-analytic at the critical point g = Qc = where the 2 solutions 
of the quadratic equation @ coalesce. There the minimum of the effective action in @ 
disappears and it becomes unbounded 0, ||, so that the saddle point solution is no longer 
valid. There is a second order phase transition at the coupling g = Qc with susceptibility 
exponent 7^°-* = | |]l], ^. This critical point is identified as the "continuum" limit of the 
random polymer theory where the number of branches, and hence the lengths, of the tree- 
graphs becomes infinite. The higher-loop contributions (molecular networks) can be found 
in [^] and their "continuum" limit is associated with an infinite number of molecules, thus 
tracing out a continuum filamentary surface |1|, |^ . 

Random polymer models can also be generated by the fermonic vector model with partition 
function 

ZF{t,g;N) = j d^pdi) e^^^-M^^f (12) 

where t/^j, i = 1, . . . , A^, are independent anticommuting nilpotent variables, ijjil) = Z^ili i^ii^i 
and the integration measure dip dip = Oili dipi dipi is defined using the usual Berezin rules 
for integrating Grassmann variables, / dipi ipi = 1, J dipi 1 = 0. The model (p!2| ) possesses 
a continuous symmetry, ipi — > J2jUijipj,ipi J2j 4'j{U''^)ji with U E GL{N,€). There is a 
further discrete symmetry under the "chiral" transformation ipi ipi, ipi ~ipi for any i, 
which is the analog of the reflection symmetry (pi —(pi of the 0{N) vector model. Together, 
these two symmetries restrict the observables of the model to those which are functions only of 
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ipip- We shall now show that the model (|T2|) possesses a random geometry interpretation and 
critical behaviour similar to that of the 0{N) vector model. A main difference with the 0{N) 
vector model is that the integration over Grassmann variables in the generating function (|T2|) 
for the polymers is a well-defined finite polynomial in the coupling constants g and t because 
nilpotency of the components ipi and ipi implies that (■?/'■?/') ^"""^ = 0. The dimension itself 
provides a cutoff on the number of terms (and polymers) in the Feynman diagram expansion 
of (|l3). Here, rather than making the partition function integration well-defined as in the 
bosonic case, the large- limit is needed to generate the full ensemble of randomly-branched 
chains. 

The Feynman diagrams for the fermion vector theory ([T2|) have propagator 



and the four-Fermi interaction vertex is 

(iJiipjipki^i") ^ = ^ {Sij6ki - Sii6jk) (14) 

These Feynman rules have the same graphical representation shown in Fig. 1 with a left- 
handed orientation for the lines (ingoing lines into a vertex for ipi components and outgoing 
lines for ipj components). Now the graphs are formed from all connected 4-point diagrams 
which preserve this orientation. The Feynman rules also associate a factor of —1 to each 
fermion loop in a Feynman graph. 

The perturbative expansion of (|12D 



Mt,0;N) ^„ n\ \2N 
is completely determined by the normalized Gaussian moments 



/f Jdij # e*'^^ \dt 

to be 



ZF{t,0;N) ' n\{N - 2n)\ \2Nt\ 

where N2 = N/2 (respectively (A^ — l)/2) when A^ is even (odd). The perturbation series is 
a finite sum which represents the same sort of random walk distribution as shown in Fig. 2 
except that it only includes polymers with up to N2 bonds. From a diagrammatic point of 
view, the alternating nature of the series arises from the minus sign associated with fermion 
loops. Term by term, this series can be made identical with the terms of same order in g in 
(|^) by the analytical continuation A^ —N/2 in (|T7p (so that Aa = A^2 in (0)). The factor of 
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2 is associated with the doubhng of degrees of freedom in the fermionic case. Thus, after the 
substitution —>■ N/2, the large expansion of the fermionic vector model is identical to 
that of the 0{N) vector model except that it is an alternating series in The coefficient 

of in the former and l/(— A^)'^ in the latter are identical. Now, however, the alternating 
nature of the fermionic vector series makes its ^ oo limit Borel summable, and as such it 
defines a better behaved statistical theory. 

To explicitly carry out the -^-expansion of the fermionic vector model, we introduce a 
scalar Hubbard-Stratonovich field (p by inserting the identity 1 = J d(p into the 

partition function integral (0) to write it as 

Zp{t,g] N)= j dip e-9^'l^^ [ d^ dij e^i-^9^/N)^^ = N\ j dip {t - igip/N)^ e^^"'^^ (18) 



When N ^ oo the integral in (p!8[) is determined by the saddle-point value of ip. Rescaling 
ip — i> ip/N , this can be found from the stationary condition for the effective action 

S{^) = -Ngp>^/2 + N\og{t - tg^) (19) 

appearing in ([T8|) which is 



tg^ - ig""^^ + ig = Q (20) 
The solution of (pO[) which is regular at (? = is 



^0 = 7^ 1-Ji-t| (21) 




Substituting (|2T|) into (|19D we get the tree-level fermionic free energy 

1 t ( n, -\ . (t 



Ff = 2 - 4^ - - j - log + V^^' " ) ^^^^ 

Note that the saddle points in ([T0|) and (0) are related by the correspondence 0^ = ^ip 
which makes more precise the analytical continuation discussed above. Notice also that the 
free energy ( P^ is related to the free energy ([TT]) of the scalar model as = 1 - 2F|°\ 
as anticipated since the large-A^ limit of the 2 models represents the same combinatorical 
problem of enumerating tree-graphs. 

The higher-loop contributions (which count the polymer networks with a given number 
of molecules) can be found by carrying out the saddle point calculation of the integral 



to higher orders. For this, we decompose the Hubbard-Stratonovich field a.s ip = ipo + ipq and 
expand the action (|T9|) in a Taylor series about the saddle-point value (^1]) in terms of the 
fiuctuation fields ipg. Using the saddle-point equation (pOD when evaluating the higher-order 
derivatives S^^\ipo), this Taylor series is found to be 

S{^) = S{ipo) -^{g + g'ipl) ^l-N± (23) 

^ n=3 ^ 
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The genus 1 free energy is then obtained from the fluctuation determinant that arises from 
Gaussian integration over the quadratic part in of (|23|), 

F^F^ = ^rM9 + g'^t) = T^r log (^9-'^ + (24) 



2N 

which also agrees with the 1-loop free energy of the 0{N) vector model 0. 

The -^-expansion of the fermionic free energy also becomes non-analytic at the critical 
point g = Qc = It exhibits the same critical behaviour as the 0^ theory above and it 

therefore lies in the same universality class as this statistical model. It is straightforward to 
carry out the double-scaling limit of the fermionic vector model in much the same way as in 
the bosonic case [0, 0, 0. This limit is associated with the continuum limit of the polymer 
network at — > oo, (7 — > in such a way that a coherent contribution from all orders of the 
perturbative and expansions is obtained. We approach the critical point Qc by defining 
a dimensionless "lattice spacing" a by at^ = 9c~ 9 and taking the continuum limit a ^ 0. 
With the rescalings mentioned above, the contribution to an arbitrary £-loop vacuum diagram 
with £ > 1 is N~^{y/a)~^N"' = N^~'^{y/ay^^^'"'. As shown in [0, the maximum number of 
vertices that a 4-point polymer diagram with i loops can have is n = 2{i — 1), so that the 
most singular behaviour of an £-loop diagram in the continuum limit a — * is (Na^^'^y~^. 
The proper continuum limit wherein a finite contribution from arbitrary genus polymers is 
obtained is thus the "double-scaling" limit where N 00 and a — in a correlated fashion so 
that the renormalized "cosmological constant" (or "linear string tension" ) A = Na?/"^ remains 
finite. The double scaling limit enables an explicit construction of the genus expansion of the 
continuum polymer theory from the vector model. 

We can now take the double scaling limit of the partition function ( |l2l) and write it as a 
loop expansion in the linear string tension A. However, as noted for the 0{N) vector model 
0, the tree and 1-loop contributions are singular in this limit. The saddle-point value (^l]) 
can be written in terms of the lattice spacing as gipo = {2/it){l — 2^/a), from which it follows 
that the genus and 1 free energies (|22|) and (0) are given by 

iVf.(»)._^_iV,og(i)+6ivVV/'-fA . iVFW-.i,ogg)+ilog(^) 

(25) 

in the continuum limit a — » 0. The A-dependent terms in (^) diverge in the double scaling 
limit and represent a non-universal behaviour of the statistical polymer system. The tree-level 
and one-loop order Feynman diagrams should therefore be subtracted in the definition of the 
double-scaling limit leading to a renormalized partition function Z/j(A,t) that only contains 
contributions from the £-loop diagrams with i >2. 

This renormalized partition function is obtained by integrating over that part of the action 
involving n > 3 vertices in the fluctuation field (fg weighted against the Gaussian form in 



7 



To pick out the finite contribution in the double-scahng hmit, we rescale the fluctuation fleld 
as yjg — > 2N^/'^ I\^/^g~^/'^ipq and note that with this rescahng we have 



in the continuum hmit a — 0. The n > A vertex terms in (^61) vanish in the double scahng 
hmit, and therefore the exact renormalized partition function in the double scaling limit is 
(up to irrelevant normalization factors) 



Ar=l,j=i , 



(27) 



The Gaussian moments in ( pTl) can be evaluated as in (g). The odd moments vanish, while 
the even moments yield a factor (3A; — 1)!!. Thus the double-scaled renormalized partition 
function admits the exact genus expansion 



Zn{A, t) = Yi-lY f^' ^V.. (28) 

' ^ ' (2^)1(512)2^ ^ ' 

The partition function (pSj) has a similar structure as that in the 0{N) vector model 
where the genus expansion is an asymptotic series with zero radius of convergence 0. In the 
fermionic case, however, the genus expansion is an alternating sum, and is therefore Borel 
summable. The convergence of the sum over genera is easily seen in the integral expres- 
sion ( p7D where the unbounded cubic term contains a factor of i which makes the overall 
integration there flnite. The Borel summability of the genus expansion is a feature unique 
to the fermionic models that does not usually occur for random geometry theories. In this 
sense, the fermionic vector model represents some novel discretized surface theory in which 
the topological expansion uniquely specifies the generating function of the statistical theory. 

The two functions 



Z±{t, g; N) = ^Zs{t, g; N)Zp{t, g; iV/2)±i (29) 

are partition functions of random polymer models where, in the case of the + sign the genera 
are restricted to be even and in the case of the - sign the genera are odd. The resulting 
statistical theory is not Borel summable, but it does represent a new sort of generating 
function for "reduced" polymer systems. In these reduced statistical models, at least part 
of the non-universal genus or 1 double-scaling behaviour is removed. The convergence 
properties of the fermionic vector models can therefore be exploited to combine them with 
other vector models and generate unusual statistical models of random polymer systems. It 
would be interesting to give the fermionic nature of these models an interpretation in terms 
of random surfaces directly. The identification e^ in the fermionic case suggests that 
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the associated random polymer theory has a complex- valued "fugacity" /i = ivr + /ig, /io G IR, 
with doubly-degenerate degrees of freedom at each vertex. Heuristically, the complex fugacity 
/i can be thought of as a rotation of the real-valued one /iq required to compensate the extra 
vertex degree of freedom that does not appear in the bosonic models. From an analytic 
point of view, the genus sum alternates relative to that of the 0{N) vector model because 
the saddle-point (^) is imaginary in the fermionic vector model ( [T^ ) so that its saddle-point 
expansion is the analytical continuation 0q = |y9o of that for the scalar model. 

The above analysis can be straightforwardly generalized to an interaction of the form 
g{ipip)^ , which will then represent a random polymer model with up to 2K-valence vertices. 
The critical behaviour is the same as that in a scalar vector model and leads to the 
same susceptibility exponent 7*^°^ = |, i.e. such a theory of random polymers is universal. To 
generate more complicated polymer models, for instance those with matter degrees of freedom 
at the vertices of the discretization |^ , one must study vector models with more complicated 
interactions. One particularly interesting aspect of the fermionic vector model is the extent 
to which these more complex models can be solved explicitly. Consider the more general 
fermionic vector model with partition function 

# # e^^^'^'^) (30) 



where V{z) is some "potential" function. The integration in (|30|) is well-defined and finite 
if e^^'-^-' has a well-defined Taylor expansion to order in the variable z. In that case, the 
partition function (^) at finite can be formally evaluated by inserting the delta function 
1 = J dz S{z — ijjip) and using the identity 



TV 



^ jdijd^ iV!(-zw;)^e™^ = N\ (-^] 6{z) (31) 



to obtain 



A^!.(|)"e-« 



z=0 



Similarly, the correlators for finite are given by 

N 



z=0 



(32) 



(33) 



The generating function and observables for these fermionic polymer models can therefore 
always be explicitly and uniquely specified. It is straightforward to check that ( p2D reduces 
to (|1^) when the potential is the four-Fermi interaction that we studied above. 

To treat the model (|30D at = 00, we use the first part of the identity (^) to write the 
partition function as 



271 



dz dw eNV{z)+^wz+N logw ^34^ 
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If the potential is a polynomial of degree m, 

m 

V{z) = Y: 9kz' (35) 

k=l 

then we can rescale z —>■ z/N and the coupling constants gk ■ simultaneously so 

that the effective action in ( P^ is NV{z) + iNwz + N\ogw. The integral ( P^ at large- is 
determined by the saddle-point value of this effective action. In the 2-dimensional complex 
space of the variables w and z, the stationary conditions are 

V'{z) +iw = Q , iz + l/w = (36) 

which can be combined into the single equation 

zV'{z) = 1 (37) 

The equation ( P7D is identical to the stationary condition for the 0{N) vector model 
defined with potential V^(0^) 0, Thus the critical behaviour of the fermionic vector model 
([30| ) is the same as that for the 0{N) vector model with the same polynomial potential (|35D . 
The genus expansion is generated by the 2-dimensional saddle-point evaluation of the integral 
([3^). The imaginary saddle-point values given by ( p6D will lead to an alternating genus 
expansion in the double-scaling limit for the fermionic theory, leading to a Borel summable 
polymer model, in contrast to the scalar case. In this case the critical point is again that 
point in coupling constant space where the function zV'{z) vanishes. For a potential of the 
form (|35D, we can adjust the coupling constants in such a way that the critical point is a 
zero of zV'{z) of order m. The leading singularity of the free energy will then be 
P, 0] which leads to the critical susceptibility exponent 7^°^ = 1 — 1/m, m = 2, 3, . . .. This 
is the multicritical series for generalized random polymer systems in dimension D > |l| 
which interpolates between the Cayley tree at m = 00 with 7*^°-' = 1 and the ordinary random 
walk we discussed earlier at m = 2 with 7*^°) = | (the A^°-component of the vector model 
free energy represents the self-avoiding random walk). In the 0{N) models, the former case 
would represent a phase of bosonic strings in target space dimension D > 1 while the latter 
case would represent a phase of pure 2-dimensional quantum gravity. In the general case, the 
potential (^) leads to discrete filamentary surfaces which have vertices of even valence up to 
2m. It would be interesting to determine what physical systems the fermionic vector models 
represent in the continuum limit. 

We see therefore that the fermionic vector model partition function is a well-defined, finite 
generating function for polymer configurations which can always (for any A^) be evaluated 
exactly. The nature of its genus expansion, i.e. that it has an alternating, Borel summable 
double scaling series, agrees with conjectures concerning the nature of the topological expan- 
sion in fermionic matrix models which have been studied in [P| as alternative random surface 
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theories to the conventional Hermitian matrix models. There it is conjectured that the genus 
expansion is an alternating series which may be Borel summable but otherwise coincides with 
the usual Painleve expansion [||. This is expected to be only true for odd polynomial po- 
tentials as it is only in that case that the matrix model possesses a chiral symmetry and 
imaginary endpoints for the support of the spectral distribution. In the case of the simpler 
fermionic vector models the partition function is always invariant under chiral transformation 
of the fermionic vector components. Furthermore, in the case of fermionic matrix models the 
correspondence with a scalar theory is more complicated - a polynomial fermion model can 
be analytically continued to a Hermitian matrix model with a generalized Penner potential 
which consists of a logarithmic plus polynomial potential fill. It should therefore represent 



a Borel summable generating function for the virtual Euler characteristics of the discretized 
moduli spaces of Riemann surfaces (rather than just the generating function for a random 
surface triangulation itself). In the vector case, the fermionic model represents the same type 
of random surface theory as the corresponding 0{N) vector field theory. Nevertheless, further 
study of the more technical features of fermionic vector models, such as an analysis of the 
Schwinger-Dyson (loop) equations and the connection with integrable hierarchies, may help 
in understanding what the corresponding structures will look like in the case of fermionic 
matrix models. In the vector case, they will generate examples of fermionic models which are 
explicitly solvable. 

We wish to thank E. Elizalde for comments on the manuscript, and S. de Queiroz and J. 
Wheater for interesting discussions. 
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